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Abstract 

In this paper, we show that the distribution of the longitudinal momentum fraction of charged hadrons dN/dz inside 
jets stemming from proton-proton collisions at -\fs = 7 TeV center of mass energy can be described by a statistical 
jet-fragmentation model. This model combines microcanonical statistics and super-statistics induced by multiplicity 
fluctuations. The resulting scale dependence of the parameters of the model turns out to be similar to what was 
observed in electron-positron annihilations in Ref. Q. 
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1. Introduction 

When calculating the spectrum of hadrons produced 
in high-energy proton-proton (pp) collisions using per- 
turbative quantum chromodynamics (pQCD) improved 
parton model calculations HI, hadron production is 
described by fragmentation functions [2, 3|. Though 
the evolution of these fragmentation functions with the 
scale Q 2 can be understood within the framework of 
pQCD El[5]|6), their actual form at a given scale Q 2 = 
s can not be deduced from QCD. In Refs. Q IS |9], it 
has been shown that fragmentation functions of charged 
hadrons, 7r-s, K-s and A-s produced in electron-positron 
(e + e~) collisions can be described by a simple statisti- 
cal physical model. This model treats hadrons formed 
in a jet as a microcanonical ensemble and obtains frag- 
mentation functions via smearing the single hadron dis- 
tribution in a jet over the distribution of charged hadron 
multiplicity measured in e + e~ annihilations. The hadron 
multiplicity distribution is an imput in this model and its 
derivation is out of the scope of this paper. Such calcu- 
lations in the microcanonical ensemble have been per- 
formed by other groups to describe hadron multiplicity 
distributions from <? + e~ to nucleus-nucleus (AA) colli- 
sions iflOl - lfT2l. 

The hadron multiplicity distribution used in Ref. Q, 
namely the Euler-Gamma distribution, belongs to the 
family of the Koba- Nielsen -Olesen (KNO) type dis- 
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tributions |fl9l 201 that describe measurements in e + e 
collisions well H2T1 l22l l23ll . 

If hadrons created in a single event are distributed 
according to the Boltzmann - Gibbs distribution, and 
the multiplicity or the temperature parameter of the 
distribution fluctuates according to the Euler-Gamma 
function, then the average hadron spectrum will be 
the Tsallis - Pareto (or q-canonical) distribution |24| - 
ll29l . Similarly, if in a single event hadrons have 
microcanonical distribution, and multiplicity fluctuates 
according to the Euler-Gamma distribution, then the 
resulting average hadron spectrum will be a micro- 
canonical generalisation of the Tsallis - Pareto (or q- 
microcanonical) distribution Q 0. It is interesting, 
that the q-microcanonical distribution can be obtained 
within the framework of non-additive thermodynamics 
too (30). 

In this paper, we point out that the fragmentation 
functions measured in pp collisions at yfs = 7 TeV 
collision energy (32] [33] can be described by the q- 
microcanonical distribution too. Furthermore, the pa- 
rameters of the distribution show similar scale depen- 
dence as seen in Refs. Il34ll35l . 

The structure of this paper is as follows: in Section[2j 
we provide a description of the microcanonical frag- 
mentation model 0. Section [3] contains the fits of the 
model to fragmentation functions and multiplicity dis- 
tributions measured in pp collision at — 7 TeV cen- 
ter of mass energy at the LHC f32l [331 - Finally, the 
conclusion is presented in Section[4] 
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2. Microcanonical Jet Fragmentation 

If the process of the creation of hadrons h\, ... , 
by the leading part on pi of a jet with multiplicity TV is 
such that the corresponding cross-section is simply pro- 
portional to the phase space available for the hadrons, 
restricted only by the energy conservation, 



\M\ 2 6 (4) 



da 



dQ. 
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then the hadrons created in the fragmentation process 
form a microcanonical ensemble. In Eq. ([T}, Q is the 
phase space of the created hadrons, p^ is the four- 
momentum, €h is the energy of the hadron hj, Pp L is the 
four-momentum of the leading parton pi, Ej el = P° pL is 
the energy of the jet, and M is the matrix amplitude 
describing the process. This way, neglecting hadron 
masses, the energy distribution of a hadron inside a jet 
with multiplicity N equals to (see Ref. J7)) 



f N (e) =A mc (l-x) 



D(N-l)-\ 



(2) 



where x = £/Ej el , e is the energy of the hadron, D 
is the effective dimensionality of the jet and A mc - 
{^ N D x )DNI(k D E° et ) follows from the normalisation 
condition 

N = j da p j dpp D -'f N (e), (3) 

with ko = f d£l p being the angular part of the mo- 
mentum space integral. Eq. |2} follows from the mi- 
crocanonical momentum space volume at fixed energy 
and multiplicity, 



Q N (E) = 



N I 
D 



• U !=1 



N \ 



N 



"5> 



7=1 > 



(4) 



N\T(DN) 

and the one-particle distribution is obtained as 

, , . Qjvr-iC^-e) 

As particles in a jet form a microcanonical ensemble, an 
entropy S j et = lnQ^i^,) and so a thermodynamic al 
temperature 
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based on the zeroth law of thermodynamics 1 30 . 31 1 can 
be associated to them. 

Microcanonical treatment of hadron production has 
also been proposed in Refs. Q - 0~8| for e + e ~, pp 
and A A reactions. The main difference between our 
approach and the ones discussed in IfTSl - Ifl8l is that 
in order to analyse the distribution of charge-averaged 
hadrons inside jets of very high energy and small jet- 
cone, we do not deal with charge conservation and ne- 
glect masses and transverse momenta of hadrons (trans- 
verse with respect to the jet direction). Thus, jet masses 

are neglected in our calculations: M 2 . ef = pr\ = 

(2 6,) 2 - (2 p,) 2 ^ 0. Consequently, the conservation 
of four-momentum is equivalent to energy conservation 
inside a one-dimensional directed jet. This way, instead 
of the jet mass Mj e , the jet energy Ej e , would control 
the distribution of hadron multiplicity. This multiplicity 
distribution we do not derive here, we rely on empirical 
fits to measurements instead. 

In Refs. l24l - ll29ll . it has been shown that special 
event-by-event fluctuation patterns of the temperature 
or of the particle multiplicity can result in power-law 
tailed average particle spectra. This applies even if in 
each event, particles are distributed according to the 
Boltzmann-Gibbs distribution. In Refs. OH - (23), it 
has been argued that an approximate Koba - Nielsen - 
Olesen (KNO) scaling of the multiplicity distribution of 
charged hadrons holds for electron-positron collisions 
(though the scaling is weakly violated by the scale evo- 
lution of the strong fine structure constant a s (Q 2 )). 

If we consider multiplicity fluctuations of the form 



p ( Ao = ^(/v-M,rv«\ 
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and microcanonical single hadron distribution inside 
each jet (cf. Eq. |2])) the multiplicity averaged hadron 
spectrum becomes 
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This result can be obtained by replacing the discrete 
sum by an integral, and using Stirling's formula n! m 



V27rn (n/e)". In terms of the integration variable £ 
X - No, only the highest power is taken into account. 

In Eq. (|8j, No is the minimal number of hadrons that 
must be produced in the fragmentation process. The 
newly introduced parameters are: q — 1 + l/(a + D + 1) 
and T = E jet /3 / [D{a + D + 1)]. The parameter q mea- 
sures the deviation of Eq. ([H} from the microcanonical 
distribution Eq. ([2}. For e + e~ data, g > 1 holds, how- 
ever, in the limit of q — > 1, the hadron distribution 



D^, (if JV/JVb -» 0) . 
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Finally, from Eqs. (Ill and ( 12 1, one can conclude that 



T = T jet [l-(D + 2)(q-l)] + O 



(,/*■) 
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I^^Ajt^d-jeyw-D-i, (9 ) 

is recovered with /V = a//3 + Nq being the mean multi- 
plicity. 

Since the multiplicity fluctuates from jet to jet, so 
does the thermodynamical temperature introduced in 
The distribution of Tj et can be obtained from 



Eq. I 

Eqs. (|6]i and Q: 



3. Analysis of Fragmentation Functions Measured 
in V* = 7 TeV Proton-Proton Collisions 

In the jet-analysis reported in Refs. 11321 [331 . very nar- 
row jet-cones of R — ^Arj 1 + A(f> 2 = 0.6 were used 
where A(f> and Arj are the azimuthal angle and pseu- 
dorapidity of the hadrons relative to that of the jet. 
(rj — - In tan 6, with 9 being the polar angle.) For such a 
jetcone, it is reasonable to make the approximation 



p{Tjet) = J^J?-(?(e-e T-ie-w- 
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with 9 = E Jet /(DT jet ), 6»o = E je ,/(DT jet0 ) and T Jet0 = 
Ejet/[D Nq-1]. The mean value of the thermodynamical 
temperature is 



z =TEvr* cosA0: 



(16) 



Furthermore, jets may be considered to be one- 
dimensional bunches of ultra-relativistic particles. This 
way, the four- momentum of the jet can be approximated 

as 
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The T parameter appearing in the multiplicity averaged 
hadron spectrum Eq. (|8]l on the other hand may be re- 
ferred to as "mean equipartition temperature". It is pro- 
portional to the average energy per particle in a jet 
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In the limit q — > 1 the mean energy per particle tends to 
the familiar result: 
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It is also worth noting that if /V » /Vo the usual equipar- 
tition formula holds for Tj et : 



P\ er = (M r coshy,M r sinhy,P7-) 



(Pt cosh rj, Pj sinh 77, P7-) 



(17) 



where Mj and P^ are the transverse energy and momen- 
tum of the jet, andy = 0.5 ln[(E jel +Pj etz )/(E j et -P jen )]. 
In the following, we will analyse jets mainly trans- 
verse to the beam direction (77 = 0), thus, we may use 
Ejet ~ Pt- Finally, the z distribution of charged hadrons 
takes the form 
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with r = T/P Tjet . 

In the canonical limit, z <s 1 Eq. Q tends to 

the Boltzmann - Gibbs distribution, and Eq. ( fl~8] > ap- 
proaches the q-canonical distribution 
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Fits of Eq. ( p~8] >, ( fT9] i and (|7]i to data on fragmentation 
functions and multiplicity distributions measured in pp 



collisions at Vs = 7 TeV l32l[33l are shown in Figs. |T|- 
|4] Figs.[T]and[3]show that the q-microcanonical (or mi- 
crocanonical Tsallis - Pareto) distribution, Eq. ( p"8j ), de- 
scribes data on dN/dz well. The low z deviation of the 
model from the data is assumably due to the low pr cut 
used in the jet analysis. Particles with transverse mo- 
mentum less than pro = 0.5 GeV/c were not taken into 
account in the jet analysis, and the downward curl of 



the measured data at low z from Eq. ( 18 i starts around 
Zo - Pro /Prjet- Similar conclusions may be drawn 
from the data-over- theory plots shown in Fig. [2] for the 
q-canonical distribution, Eq. ( 19 1, except that this dis- 
tribution describes data for z ~ 0.2 only. 

The evolution of the fitted q and T parameters of 
Eqs. (JT8]l and (fT9~ll with the transverse momentum of 



the jet are shown in Figs. [5] - [8] The "mean equipar- 
tition temperature" parameter scaled by the transverse 
momentum of the jet T* = T /Prjet shows power-law 
dependence on Prjet, 

T* = (PrjetfOof . (2°) 
while for the parameter q, both a power-law, 



q = (Prjet/Qof 
and a double-logarithmic ansatz, 



q = 1 + y ulnln(p 7> ,/go) 



(21) 



(22) 



fit. The 2o an d I 1 parameters of the q-canonical and the 
q-microcanonical distributions approximately coincide 
within errors. 

For q, Eq. ( 22 1 was successfully used in Ref. l35l to 



adjust the Q 2 evolution of a fragmentation function of 
the form of Eq. ([19]) to that of an AKK type one Q. 
In pp collisions, Q = Prjet seems to be a good choice. 
The q and T parameters show similar scale dependence 
for fragmentation functions of protons, K°-s, n°-s, A-s 
and charge-averaged hadrons produced in e + e~ annihi- 
lations as well as for the transverse momentum spec- 
tra of charged hadrons stemming from pp collisions 

HE EM. 

Fig.|4]shows that Eq. <|Vj describes data on multiplic- 
ity distributions well, except for N < 3, where measured 
data are higher than what the Euler-Gamma distribution 
predicts. This effect is perhaps due to the small cone- 
size. In a jet with one or two very energetic particles, 
the others have very small energies, and thus may fly 
out of the jetcone. This way, the number of jets with 
only a few particles increases in this type of jet analysis. 
As a consequence, the a parameter of Eq. which 
is greatly influenced by the number of low multiplicity 



jets, can not be determined reliably. Thus, it is not so 
disconcerting that the q parameter predicted from mul- 
tiplicity fits takes lower values (of around q = 1.1) than 
that obtained from fits to dN/ dz data. 

It is important to note that in Ref. [33], jets were 
reconstructed from charged particles only, while in 
Ref. 1321 . calorimetric measurements were used in the 
jet reconstruction, so both neutral and charged particles 
were included in the analysis. For this reason, we fitted 



Eqs. (20 1 - (22 1 to q and T values obtained for jets with 
high transverse momenta only (25 GeV/c < Prjet ^ 500 
GeV/c from l32l ). Nevertheless, q and T values for jets 
with low transverse momenta (4 GeV/c < Prjet - 40 
GeV/c obtained from [33 1) show a tendency similar to 
that of the high Prjet jets. The T* parameters of the low 
Prjet dataset are approximately 10% higher, than what 



the fit of Eq. ( 20 1 to the high Prjet dataset predicts (see 



Pig. [8j». From Eq. ( 12 1, it can be seen that T* is propor- 
tional to the inverse of the multiplicity. This way, the 
ratio of T*-s obtained in the two different analysis' is 
proportional to the ratio of the total multiplicity to that 
of charged particles. If as an estimate, we used the ra- 
tio of charged to neutral pions, we would get a factor of 
3/2 for the ratio of T*-s obtained from the two different 
analysis'. This value is somewhat higher than what can 
be seen in Fig. [8] 

Eq. ( 18 1 describes both datasets. For the high Prjet 
dataset, the power of the 1 — z factor in the numerator 
takes the value D(No - 1) - 1 = +1, while this quantity 
decreases from 1 to - 1 for the low Prjet dataset as Prjet 
decreases from 40 GeV/c to 4 GeV/c. 

4. Conclusions 

This paper shows that the statistical jet-fragmentation 
model 1 7 , 8 1 describes the dN/dz distribution of hadrons 
in jets created in proton-proton reactions at V* = 7 TeV 
center of mass energy Il32l [33l . This model combines 
microcanonical statistics (which has also been used 
in the description of different hadronic observables in 
high-energy phenomena in Refs. |[T0ll - lfT71l ) with super- 
statistics [24|-|29| stemming from multiplicity fluctua- 
tions emerging in proton-proton as well as in electron- 
positron and nucleus-nucleus collisions 1 19 1-[23 1. 

It turns out that the parameters of the dN/dz dis- 
tribution of charged hadrons in jets in proton-proton 
collisions (Sect. [3]l show similar scale dependence as 
the parameters of fragmentation functions in electron- 
positron annihilations Q |8] [9j and of transverse mo- 
mentum spectra of charged hadrons stemming from 
proton-proton collisions |34|. These scale evolutions 
are consistent with the DGLAP equations ||35l . 
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Figure 1 : Measured distributions of the longitudinal momentum frac- 
tion z of hadrons inside jets with various transverse momenta (data 
of jets with P Tje , = [4 - 6] , . . . , [24 - 40] GeV/c and with P Tjel = 
[25 - 40] , ... , [400 - 500] GeV/c are published in Ref. (33] and in 
Ref. 1321 respectively) and fitted 1 dimensional q-microcanonical dis- 
tributions (Eq. (18) with D = 1 , and No = 3 for the high Prjei dataset 
and Nq = 1, 1, 1, 2,3 (from bottom to top) for the low Pjjet dataset). 




Figure 3: Ratios of measured dN/dz distributions and fitted 1 dimen- 
sional q-microcanonical distributions (Eq. Q8) with D = 1. For the 
values of No, see the caption of Fig.[TJ for jets with various transverse 
momenta (data of graphs are published in Refs. I33II32I ). 




N/n 



Figure 2: Ratios of measured dN/dz distributions and fitted 1 di- 
mensional q-canonical distributions (Eq. j!9) with D = 1) for jets 
with various transverse momenta (data of graphs are published in 
Refs. f33l[32l ). 



Figure 4: Measured multiplicity distributions of charged hadrons in- 
side jets with various transverse momenta and rapidity and fitted 
Euler-Gamma distributions (Eq. ^J). Data of graphs are published 
in Refs. 13311321 . 
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q = 1 + {p / Q/ 

Q = 28398 ± 17682 GeV/c 
H = 0.15 ± 0.02 




q = 1 + ji log log{p T]et /Q Q ) 
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H = 0.32 ± 0.02 
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Figure 5: Fitted values of the q parameter in Eq. l |19( with D = 1 to 
measured dN/dz distributions shown in Fig.[T] 



Figure 7: Fitted values of the q parameter in Eq. j!8) with D = 1 
(for the values of No, see the caption of Fig.^ to measured dN/dz 
distributions shown in Fig.|T| 
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Figure 6: Fitted values of the T* parameter in Eq. j!9) with D = 1 to 
measured dN/dz distributions shown in Fig.[T] 



Figure 8: Fitted values of the T* parameter in Eq. jl S[ with D = 1 
(for the values of No, see the caption of Fig. [TJ to measured dN/dz 
distributions shown in Fig.^ 
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Finally, it is pointed out that the dN/dz distributions 
obtained from two different jet analysis' and different 
kinematical ranges f32l [33l both can be described by 
the microcanonical jet-fragmentation model. 
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